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Abstract 



We report new results on nearly conformal gauge theories with fermions in the fundamental representation of the SU(3) color gauge 
group as the number of fermion flavors is varied in the Nf = 4-16 range. To unambiguously identify the chirally broken phase 
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below the conformal window we apply a comprehensive lattice tool set in finite volumes which includes the test of Goldstone 
pion dynamics, the spectrum of the fermion Dirac operator, and eigenvalue distributions of random matrix theory. We also discuss 
the theory inside the conformal window and present our first results on the running of the renormalized gauge coupling and the 
renormalization group beta function. The importance of understanding finite volume zero momentum gauge field dynamics inside 
the conformal window is illustrated. Staggered lattice fermions are used throughout the calculations. 
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1. Introduction 

The Large Hadron Collider will probe the mechanism of 
electroweak symmetry breaking. It is an intriguing possibil- 
ity that new physics beyond the Standard Model might take 
the form of some new strongly-interacting gauge theory. In 
one scenario, the Higgs sector of the electroweak theory is re- 
placed by a so-called technicolor theory, whose dynamics pro- 
vides the required spontaneous symmetry breaking Q] [2] [3]. 
These models avoid the fine-tuning problem and may lead to a 
heavy composite Higgs particle on the TeV scale. Although 
attractive, the challenge is to extend a technicolor theory to 
include fermion mass generation, while satisfying the various 
constraints of electroweak phenomenology. This idea has lately 
been revived by new explorations of the multi-dimensional the- 
ory space of nearly conformal gauge theories [5] [6] [7) . The 
terminology of technicolor in this report will refer in a generic 
sense to these investigations. Exploring the new technicolor 
ideas has to be based on nonperturbative studies which are only 
becoming feasible now with the advent of new lattice technolo- 
gies. 

Model building of a strongly interacting electroweak sector 
requires the knowledge of the phase diagram of nearly con- 
formal gauge theories as the number of colors N c , number of 
fermion flavors Nf, and the fermion representation R of the 
technicolor group are varied in theory space. For fixed N c and R 
the theory is in the chirally broken phase for low Nf and asymp- 
totic freedom is maintained with a negative /3 function. On the 
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other hand, if Nf is large enough, the ft function is positive for 
all couplings, and the theory is trivial. If the regulator cut-off is 
removed, we are left with a free non-interacting continuum the- 
ory. There is some range of Nf for which the f3 function might 
have a non-trivial zero, an infrared fixed point, where the the- 
ory is in fact conformal IfTOlfTTII . This method has been refined 
by estimating the critical value of Nf, above which spontaneous 
chiral symmetry breaking no longer occurs Ifl2l IT3l [141 . 

Interesting models require the theory to be very close to, but 
below, the conformal window, with a running coupling which is 
almost constant over a large energy range lfl5l [TBI [T71 [T8l [T9l . 
The nonperturbative knowledge of the critical NJ" separating 
the two phases is essential and this has generated much interest 
and many new lattice studies [20, 21 , 22, 23, 24, 25, 26 27ll28l 
|29l|3aiSl32l[^[Ml[35l[36ll37l!38, 39 40] ED @! [43] HUES 



Our goal to unambiguously identify the chirally broken phase 
below the conformal window requires the application and test- 
ing of a comprehensive lattice tool set in finite volumes which 
includes the test of Goldstone pion dynamics, the spectrum of 
the fermion Dirac operator, and eigenvalue distributions of Ran- 
dom Matrix Theory (RMT). Inside the conformal window we 
investigate the running coupling and the f3 function. We report 
new results at Nf = 4, 8, 9, 12, 16 for fermions in the fundamen- 
tal representation of the SU(3) technicolor gauge group. We 
find Nf = 4, 8, 9 to be in the chirally broken phase and Nf = 16 
is consistent with the expected location inside the conformal 
window. To resolve the Nf = 12 phase from our simulations 
will require further analysis. 
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2. Chiral symmetry breaking below the conformal window 

We will identify in lattice simulations the chirally broken 
phases with Nj = 4,8,9 flavors of staggered fermions in the 
fundamental SU(3) color representation using finite volume 
analysis. The staggered fermions are deployed with a special 
6-step exponential (stout) smearing procedure B71 in the lattice 
action to reduce well-known cutoff effects with taste breaking 
in the Goldstone spectrum. The presence of taste breaking re- 
quires a brief explanation of how staggered chiral perturbation 
theory is applied in our analysis. The important work of Lee, 
Sharpe, Aubin and Bernard (48] |49j [50j is closely followed in 
the discussion. 

2.1. Staggered chiral perturbation theory 



Starting with the Nf — 4 example [48], the spontaneous 
breakdown of S U (4)l x S U(4)r to vector S U (4) gives rise to 
15 Goldstone modes, described by fields <\>,. These can be orga- 
nized into an S U{4) matrix 



^) = «p(j__), = £ 



(1) 



where F is the Goldstone decay constant in the chiral limit and 
the normalization T a = jf^y,^} is used for the flavor 

generators. The leading order chiral Lagrangian is given by 



L 



<4) 



Tr(d„Z<9„E t ) - -Bm„ F 2 Tr(£ + S f ) , 



(2) 



with the fundamental parameters F and B measured on the tech- 
nicolor scale Atc which replaced Aqcd in the new theory. Ex- 
panding the chiral Lagrangian in powers of <p one finds 15 de- 
generate pions with masses given by 



Ml = 2Bm q [l + 0(m q /A TC )] 



(3) 



The leading order term is the tree-level result, while the cor- 
rections come from loop diagrams and from higher order terms 
in the chiral Lagrangian. The addition of a 2 -£| 6) breaks chiral 
symmetry and lifts the degeneracy of the Goldstone pions. Cor- 
rection terms are added to Eq. (|3]l which becomes 

M 2 = C(T a )-a 2 A\ c +2Bm q [l + 0(m ? /A TC ) + 0(a 2 A 2 c )] (4) 

where the representation dependent C(T a ) is a constant of order 
unity. Contributions proportional to a 2 are due to x!^\ and lead 
to massive Goldstone pions even in the m q — > chiral limit. 
The only exception is the pion with flavor £5 which remains 
massless because the U(1)a symmetry is protected. 

Lee and Sharpe observe that the part of X.f } without deriva- 
tives, defining the potential r vi , is invariant under flavor S (9(4) 
transformations and gives rise to the a 2 term in M 2 . Terms in 
£,f^ involving derivatives break S (9(4) further down to the lat- 
tice symmetry group and give rise to non-leading terms propor- 



tional to a 2 m and a 4 . The taste breaking potential is given by 

-<Vf = CiTr&E&tf) 

+ C 2 |[Tr(E 2 )-Tr(^ 5 2) + /i.c.] 
+ C 3 |^[Tr(^ v S)+/z. c .] 

V 

+ C 4 i2[ Tr ^5^5v£)+/ 1 .C.] 

V 

+ C5-2J] [TrtfyZfytf) - Tr&aEftytf)] 

V 

+ C 6 ^Tr(^ v/J E t ). (5) 

//<v 

The six unknown coefficients C, are all of size A^, c . 

In the continuum, the pions form a 15-plet of flavor S U{4), 
and are degenerate. On the lattice, states are classified by the 
symmetries of the transfer matrix and the Goldstone pions fall 
into 7 irreducible representations: four 3-dimensional represen- 
tations with flavors £5, £y and £,4, and three 1-dimensional 
representations with flavors £4, £45 and £5. 

Close to both the chiral and continuum limits, the pion 
masses are given by 

M„(T a ) 2 = 2Bm q + a 2 A(T a ) + 0(a 2 m q ) + 0(a 4 ) , (6) 

with A(T a ) ~ A 4 C arising from r V / | 6) . Since r V / j 6) respects flavor 
SO{4), the 15 Goldstone pions fall into 5(9(4) representations: 



A(ft) = 0, 



(7) 



A(&) = ^(d +C2 + C3 + 3C4 + C5+3C6), (8) 



A(^s) = ^(Ci +C 2 + 3C 3 +C 4 -C 5 +3C 6 ), 



F 2 



(2C 3 + 2C 4 + 4C 6 ) . 



(9) 



(10) 



In the chiral limit at finite lattice spacing the lattice irreducible 
representations with flavors and £4 are degenerate, those with 
flavors £5 and £45, and those with flavors and £,4 are degen- 
erate as well. No predictions can be made for the ordering or 
splittings of the mass shifts. We also cannot predict the sign 
of the shifts, although our simulations indicate that they are 
all positive with the exponentially smeared staggered action we 
use. This makes the existence of an Aoki phase ll48l unlikely. 

The method of |48 1 has been generalized in a nontrivial way 
to the Nf > 4 case 11491 l50l which we adopted in our calcula- 
tions with help from Bernard and Sharpe. The procedure cannot 
be reviewed here but it will be used in the interpretation of our 
Nf = 8 simulations. 

2.2. Finite volume analysis in the p-regime 

Three different regimes can be selected in simulations to 
identify the chirally broken phase from finite volume spectra 
and correlators. For a lattice size L 3 S x L, in euclidean space and 
in the limit L, » L s , the conditions F„L S > 1 and M„L S > 1 
select the the p-regime, in analogy with low momentum count- 
ing EUEa. 
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For arbitrary Nf, in the continuum and in infinite volume, 
the one-loop chiral corrections to M„ and F„ of the degenerate 
Goldstone pions are given by 



Mi = M- 



1 



M 



8n 2 N f F 2 



4-) 

\MJ 



1 + 



NfM 2 



M-) 

\M) 



(11) 



(12) 



where M 2 = 2B 



l6n 2 F 2 

and F, B, A3, A4 are four fundamental pa 



rameters of the chiral Lagrangian, and the small quark mass m q 
explicitly breaks the symmetry [53 1. The chiral parameters F, B 
appear in the leading part of the Lagrangian in Eq. Q, while 
A3, A4 enter in next order. There is the well-known GMOR re- 
lation 'Lcond = BF 2 in the m q — > limit for the chiral condensate 
per unit flavor Il54ll . It is important to note that the one-loop cor- 
rection to the pion coupling constant F„ is enhanced by a factor 
N 2 compared to M 2 . The chiral expansion for large Nf will 
break down for F„ much faster for a given M n IF n ratio. 

The finite volume corrections to M n and F„ are given in the 
p-regime by 



M„(L s ,ij) = M n 



F n (L s , rf) = F„ 



1 + 



1 - 



1 



M 1 



2N f l6n 2 F 2 



N f M 2 _ 



giU, rf) 



(13) 



(14) 



2 16;T 2 F 2 

where ~g\ {A, rf) describes the finite volume corrections with A = 
M ■ L s and aspect ratio 77 = L,/L s . The form of ~gi(A, rf) is a 
complicated infinite sum which contains Bessel functions and 
requires numerical evaluation ll52l . Eqs. ( TTp4 i provide the 
foundation of the p-regime fits in our simulations. 

2.3. 5-regime and e-regime 

At fixed L s and in cylindrical geometry L t /L s » 1, a 
crossover occurs from the p-regime to the 5-regime when m q — > 
0, as shown in Fig.[T] The dynamics is dominated by the rotator 
states of the chiral condensate in this limit lf55l which is charac- 
terized by the conditions FL S > 1 and ML S <K 1 . The densely 
spaced rotator spectrum scales with gaps of the order ~ \IF 2 L? S , 
and at m q — the chiral symmetry is apparently restored. How- 
ever, the rotator spectrum, even at m q — in the finite volume, 
will signal that the infinite system is in the chirally broken phase 
for the particular parameter set of the Lagrangian. This is of- 
ten misunderstood in the interpretation of lattice simulations. 
Measuring finite energy levels with pion quantum numbers at 
fixed L s in the m q — > limit is not a signal for chiral symmetry 
restoration of the infinite system ll36l . 

If L, ~ L s under the conditions FL S > 1 and ML S <K 1 , the 
system will be driven into the e-regime which can be viewed 
as the high temperature limit of the (5-regime quantum rotator. 
Although the 5-regime and e-regime have an overlapping re- 
gion, there is an important difference in their dynamics. In the 
5-regime of the quantum rotator, the zero spatial momentum 
of the pion field U(x) dominates with time-dependent quantum 
dynamics. The e-regime is dominated by the four-dimensional 
zero momentum mode of the chiral Lagrangian. 



chiral p-regime 




Figure 1: Schematic plot of the regions in which the three low energy chi- 
ral expansions are valid. The vertical axis shows the finite temperature scale 
(euclidean time in the path integral) which probes the rotator dynamics of the 
(5-regime and the e-regime. The first two low lying rotator levels are also shown 
on the vertical axis for the simple case of Nf = 2. The fourfold degenerate 
lowest rotator excitation at m q = will split into an isotriplet state (lowest en- 
ergy level), which evolves into the p-regime pion as m q increases, and into an 
isosinglet state representing a multi-pion state in the p-regime. Higher rotator 
excitations have similar interpretations. 
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Figure 2: The crossover from the p-regime to the (5-regime is shown for the n 
and 7T(5 states at Nf = 4. 

We report simulation results of all three regimes in the chi- 
rally broken phase of the technicolor models we investigate. 
The analysis of the three regimes complement each other and 
provide cross-checks for the correct identification of the phases. 
First, we will probe Eqs.( TT][l4 i in the p-regime, and follow 
with the study of Dirac spectra and RMT eigenvalue distribu- 
tions in the e-regime. The spectrum in the (5-regime is used as a 
signal to monitor p-regime spectra as m q decreases. Fig.|2]is an 
illustrative example for this crossover in our simulations. 



3. Simulations results in the p-regime 

The tree level improved Symanzik gauge action was used in 
our simulations. The link variables in the staggered fermion 
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Figure 3: The first two rows of the composite figure show Nf = 4 simulation results in the p-regime. The first row de picts th e collapsing pion spectrum and the 
techni-rho as the continuum limit is approached. The second row shows the chiral fits to M^/m q and F n based on Eqs. jl l|l4 . The range m q = 0.008 - 0.025 is 



used in the fitting procedure. The approximately linear behavior of the chiral condensate (l/nfi) is also shown in the second row. The third and fourth rows summarize 
the simulation results for Nf = 8. The third row shows the collapsing pion spectrum and the techni-rho as the continuum limit is approached. The chiral fit to 
M^lm q is shown based on Eq. 1111 only since the F n data points are outside the convergence range of the chiral expansion. The range m q = 0.015 - 0.035 is used 
in the fitting procedure. The fifth row illustrates our first simulation results for Nf = 9. It shows the split pion spectrum, chiral fit to M%/m q and the F n data points 
are outside the convergence range of the chiral expansion. The range m q = 0.02 - 0.04 is used in the fitting procedure. 
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matrix were exponentially smeared with six stout steps at Nf = 
4, 8 and four stout steps at Nf = 9. The RHMC algorithm was 
deployed in all runs but rooting of the fermion determinant only 
affected the Nf = 9 simulations. The results shown in Fig.[3]are 
from the p-regime of the chirally broken phase with the condi- 
tions M„ ■ L s » 1 and F„ ■ L ~ 1 when the chiral condensate 
begins to follow the expected behavior of infinite volume chi- 
ral perturbation theory from Eqs. ( 1 1|12 1 with calculable finite 
volume corrections from Eqs. 
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The Nf = 4 simulations work in the p-regime as expected. 
The pion spectrum is clearly separated from the technicolor 
scale of the p-meson whose quadratic fit is just to guide the eye. 
Moving towards the continuum limit with increasing B = 6/g 2 , 
we see the split pion spectrum collapsing onto the true Gold- 
stone pion. The true Goldstone pion and two additional split 
pion states are shown. A is the measure of the small quadratic 
pion mass splittings in lattice units. Their origin was discussed 
in Section 2 in Eqs. (TfTOi. The spectrum is parallel and the 
gaps appear to be equally spaced consistent with the earlier ob- 
servation in QCD where the C\ term seems to dominate taste 
breaking accounting for the equally spaced pion levels [48|. 



The simultaneous chiral fit of M*/m q and F„ based on Eqs. ( 1 1 



14 1 works when the chiral loop term corrects the tree level value 
of Ml /m q = IB. This is a chirally broken phase and the picture 



holds in the m„ 



limit. The fit to determine the N 



/ 



chiral condensate for m q — is shown in the second row on the 
right. It sets the scale of electroweak symmetry breaking in the 
Higgs mechanism. 

As we move to the Nf — 8 p-regime simulations summarized 
in the third and forth rows of Fig.|3]we observe the weakening 
of the chiral condensate and increased difficulties in passing the 
chiral tests. The pion spectrum is still clearly separated from 
the technicolor scale of the p-meson. Moving towards the con- 
tinuum limit with increasing B = 6/g 2 , we see the split pion 
spectrum collapsing toward the true Goldstone pion with a new 
distinguished feature. The true Goldstone pion and two addi- 
tional split pion states are shown with different slopes as m q 
increases. Towards m q = the pion spectrum is collapsed at 
fixed gauge coupling, indicating that the effects of leading or- 
der taste breaking operators, the generalization of those from 
Nf = 4 to Nf = 8 as discussed in Section 2, are smaller than at 
Nf = 4 in the explored coupling constant range. This is some- 
what unexpected and unexplained. Next to leading order taste 
breaking operators are responsible for the spread of the slopes 
and they seem to dominate. They were identified in Eq. (|6]l 
as the last two terms. It is reassuring to see that this structure 
is collapsing as we move toward the continuum limit. We an- 
alyzed this pattern within staggered perturbation theory in its 
generalized form beyond four flavors [49 50 1. The simultane- 
ous chiral fit of M%/m q and F„ based on Eqs. ( ll|l4 l cannot 
be done at Nf — 8 within the reach of the largest lattice sizes 
we deploy since the value of aF is too small even at L=24 for 
coupling constants where taste breaking drops to an acceptable 



level. The chiral fit to M 2 /m q is shown based on Eq. 



ccepi 



only 



since the F„ data points are outside the convergence range of 
the chiral expansion. We would need much bigger lattices to 
drop further down in the p-regime with m q to the region where 



the simultaneous fit could be made. It is also important to note 
that the chiral condensate is very small in the m q — > limit in 
the region where taste breaking is not large. This is shown in 
row four of Fig.[3]on the right side. 

The Nf = 9 p-regime simulations are summarized in the fifth 
row of Fig.|3]where we observe the continued weakening of the 
chiral condensate and the increased difficulties in passing the 
chiral tests. The pion spectrum is still clearly separated from 
the technicolor scale of the p-meson. Moving towards the con- 
tinuum limit to see the split pion spectrum collapsing toward 
the true Goldstone pion is increasingly difficult. The true Gold- 
stone pion and two additional split pion states are shown again 
with different slopes as m q increases. Forcing the collapse of 
the split pion spectrum will require larger lattices with smaller 
gauge couplings. The trends and the underlying explanation is 
very similar to the Nf = 8 case. The chiral fit to M 2 lm q is 



shown based on Eq. (Hi only since the F„ data points are out- 
side the convergence range of the chiral expansion. 

In summary, we have shown that according to p-regime tests 
the Nf = 4,8,9 systems are all in the chirally broken phase 
close to the continuum limit. Currently we are investigating 
Nf = 9, 10, 1 1, 12 on larger lattices to determine the lower edge 
of the conformal window. Lessons from the Dirac spectra and 
RMT to complement p-regime tests are discussed in the next 
section including comments about the controversial Nf = 12 
case. 

4. Epsilon regime, Dirac spectrum and RMT 

If the bare parameters of a gauge theory are tuned to the s- 
regime in the chirally broken phase, the low-lying Dirac spec- 
trum follows the predictions of random matrix theory. The cor- 
responding random matrix model is only sensitive to the pat- 
tern of chiral symmetry breaking, the topological charge and the 
rescaled fermion mass once the eigenvalues are also rescaled by 
the same factor £ eo „</ V. This idea has been confirmed in various 
settings both in quenched and fully dynamical simulations. The 
same method is applied here to nearly conformal gauge models. 

The connection between the eigenvalues A of the Dirac oper- 
ator and chiral symmetry breaking is given in the Banks-Casher 
relation l56l. 



Wand = ~m = lim lim lim 

i->0 m->0 V^oc 



V ' 



where "L con d designates the quark condensate normalized to a 
single flavor. To generate a non-zero density p(0), the small- 
est eigenvalues must become densely packed as the volume in- 
creases, with an eigenvalue spacing AA ~ l/p(0) = n/(I. con dV). 
This allows a crude estimate of the quark condensate Z com /. One 
can do better by exploring the e-regime: If chiral symmetry is 
spontaneously broken, tune the volume and quark mass such 
that j- <K L «: -g-, so that the pion is much lighter than the 
physical value, and finite-volume effects are dominant as we 
discussed in Section 2. The chiral Lagrangian of Eq. Q is 
dominated by the zero-momentum mode from the mass term 
and all kinetic terms are suppressed. In this limit, the distribu- 
tions of the lowest eigenvalues are identical to those of random 
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Figure 4: From simulations at Nf = 4 the first row shows the approach to quartet degeneracy of the spectrum as ji increases. The second row shows the integrated 
distribution of the two lowest quartets averaged. The solid line compares this procedure to RMT with Nf = 4. 



matrix theory, a theory of large matrices obeying certain sym- 
metries ll57l l58l l59l . To connect with RMT, the eigenvalues 
and quark mass are rescaled as z = XL conl jV and = m q I. coni iV, 
and the eigenvalue distributions also depend on the topological 
charge v and the number of quark flavors Nf. RMT is a very 
useful tool to calculate analytically all of the eigenvalue dis- 
tributions. The eigenvalue distributions in various topological 
sectors are measured via lattice simulations, and via compar- 
ison with RMT, the value of the condensate can be ex- 
tracted. 

After we generate large thermalized ensembles, we calculate 
the lowest twenty eigenvalues of the Dirac operator using the 
PRIMME package |60|. In the continuum limit, the staggered 
eigenvalues form degenerate quartets, with restored taste sym- 
metry. The first row of Fig. [4] shows the change in the eigen- 
value structure for Nf = 4 as the coupling constant is varied. At 
P — 3.6 grouping into quartets is not seen, the pions are notice- 
ably split, and staggered perturbation theory is just beginning to 
kick in. At f3 — 3.8 doublet pairing appears and at /3 = 4.0 the 
quartets are nearly degenerate. The Dirac spectrum is collapsed 
as required by the Banks-Casher relation. In the second row we 
show the integrated distributions of the two lowest eigenvalue 
quartet averages, 

f p k (A')dA', k =1,2 (15) 
Jo 

which is only justified close to quartet degeneracy. All low 
eigenvalues are selected with zero topology. To compare with 
RMT, we vary // = m q I, com iV until we satisfy 

vU)sim (Zl)RMT 

= . (16) 

m fi 

where (/li) s i m is the lowest quartet average from simulations and 



the RMT average (z)rmt depends implicitly on yU and Nf. With 
this optimal value of /i, we can predict the shapes of pk{A) and 
their integrated distributions, and compare to the simulations. 
The agreement with the two lowest integrated RMT eigenvalue 
shapes is excellent for the larger {3 values. 




0.01 0.02 0.03 0.04 0.05 

a a 



Figure 5: The solid lines compare the integrated distribution of the two lowest 
quartet averages to RMT predictions with Nf = 8. 

The main qualitative features of the RMT spectrum are very 
similar in our Nf = 8 simulations as shown in Fig. [5] One 
marked quantitative difference is a noticeable slowdown in re- 
sponse to change in the coupling constant. As j3 grows the 
recovery of the quartet degeneracy is considerably delayed in 
comparison with the onset of p-regime Goldstone dynamics. 
Overall, for the Nf — 4, 8 models we find consistency between 
the p-regime analysis and the RMT tests. Earlier, using Asqtad 
fermions at a particular (3 value, we found agreement with RMT 
even at Nf - 12 which indicated a chirally broken phase [20 1. 
Strong taste breaking with Asqtad fermion leaves the quartet 
averaging in question and the bulk pronounced crossover of the 
Asqtad action as /3 grows is also an issue. Currently we are 
investigating the RMT picture for Nf = 9, 10, 11, 12 with our 
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much improved action with four and six stout steps. This action 
shows no artifact transitions and handles taste breaking much 
more effectively. Firm conclusions on the Nj ■ — 12 model will 
require continued investigations. 

5. Inside the conformal window 

We start our investigation and simulations of the conformal 
window at Nf = 16 which is the most accessible for analytic 
methods. We are particularly interested in the qualitative behav- 
ior of the finite volume spectrum of the model and the running 
coupling with its associated beta function which is expected to 
have a weak coupling fixed point around g* 2 m 0.5, as estimated 
from the scheme independent two-loop beta function |61 1. 



For S U(3), T x = A 3 /2 and T 2 = A & /2. With N f flavors of mass- 
less fermion fields the effective potential of the constant mode 
is given by 

^eff(C") = J] V(C b [p.f - $ ]) -N f J] V(C»$ + TTk), (17) 

i>j i 

with k = for periodic, or k = (1,1,1), for anti-periodic 
boundary conditions on the fermion fields. The function V(C) is 
the one-loop effective potential for Nf = and the weight vec- 
tors are determined by the eigenvalues of the abelian gener- 
ators. For SU(3) ju (1) = (l,l,-2)/Vl2and// 2) = 1(1,-1,0). 
The correct quantum vacuum is found at the minimum of 
this effective potential which is dramatically changed by the 
fermion loop contributions. The Polyakov loop observables re- 



5.7. Conformal dynamics infinite volume 

A distinguished feature of the Nf = 16 conformal model is 
how the renormalized coupling g 2 (L) runs with L, the linear 
size of the spatial volume in a Hamiltonian or Transfer Matrix 
description. On very small scales the running coupling g 2 (L) 
grows with L as in any other asymptotically free theory. How- 
ever, g 2 (L) will not grow large, and in the L — > oo limit it will 
converge to the fixed point g* 2 which is rather weak, within the 
reach of perturbation theory. There is nontrivial small volume 
dynamics which is illustrated first in the pure gauge sector. 

At small g 2 , without fermions, the zero momentum compo- 
nents of the gauge field are known to dominate the dynam- 
ics (62] [63] |64). With SU(3) gauge group, there are twenty 
seven degenerate vacuum states, separated by energy barri- 
ers which are generated by the integrated effects of the non- 
zero momentum components of the gauge field in the Born- 
Oppenheimer approximation. The lowest energy excitations of 
the gauge field Hamiltonian scale as ~ g 2 ^ 3 (L)/L evolving into 
glueball states and becoming independent of the volume as the 
coupling constant grows with L. Nontrivial dynamics evolves 
through three stages as L grows. In the first regime, in very 
small boxes, tunneling is suppressed between vacua which re- 
main isolated. In the second regime, for larger L, tunneling 
sets in and electric flux states will not be exponentially sup- 
pressed. Both regimes represent small worlds with zero mo- 
mentum spectra separated from higher momentum modes of 
the theory with energies on the scale of 2n/L. At large enough 
L the gauge dynamics overcomes the energy barrier, and wave 
functions spread over the vacuum valley. This third regime is 
the crossover to confinement where the electric fluxes collapse 
into thin string states wrapping around the box. 

It is likely that a conformal theory with a weak coupling fixed 
point at Nf = 16 will have only the first two regimes which 
are common with QCD. Now the calculations have to include 
fermion loops [65, 66]. The vacuum structure in small enough 
volumes, for which the wave functional is sufficiently localized 
around the vacuum configuration, remains calculable by adding 
in one loop order the quantum effects of the fermion field fluctu- 
ations. The spatially constant abelian gauge fields parametriz- 
ing the vacuum valley are given by A,(x) = T a C/L where T a 
are the (N-l) generators for the Cartan subalgebra of SU(N). 




Figure 6: Polyakov loop distributions, blue in the time-like and red in the space- 
like directions, from our Nf = 16 simulation with 16 4 volume at /? = 18 with 
tree level Symanzik improve gauge action and staggered fermions with six stout 
steps. The fermion boundary condition is anti-periodic in the time direction and 
periodic in the spatial directions. 

main center elements at the new vacuum configurations with 
complex values 

Pj = i tr (exp(/C*r„)) = - J] ex P 0/4 n) C*) = exp&rUj/N), 

n 

(18) 

for S U(N). This implies that f/^C b = 2nl/ N (mod In), inde- 
pendent of n, and V e k ff = -N f NV(2nl/N + nk). In the case of 
anti-periodic boundary conditions, k = (1, 1, 1), this is mini- 
mal only when 1 = (mod 2n). The quantum vacuum in this 
case is the naive one, A — (Pj = 1). In the case of peri- 
odic boundary conditions, k = 0, the vacua have 1 + 0, so 
that Pj correspond to non-trivial center elements. For SU(3), 
there are now 8 degenerate vacua characterized by eight dif- 
ferent Polyakov loops, Pj = exp(+27n'/3). Since they are re- 
lated by coordinate reflections, in a small volume parity (P) and 
charge conjugation (C) are spontaneously broken, although CP 
is still a good symmetry 11651 . As shown in Fig. [6] our simu- 
lations in the Nf = 16 model near the fixed point g* 2 confirm 
this picture. In the weak coupling phase of the conformal win- 
dow the time-like Polyakov loop takes the real root, while the 
space-like Polyakov loops always take the two other complex 
values, as expected on the basis of the above picture. Next we 
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will describe our method to probe the running coupling inside 
the conformal window. It is a pilot study for more compre- 
hensive investigations of weak and strong coupling conformal 
dynamics. 

5.2. Running coupling and beta function 

Consider Wilson loops W(R, T, L), where R and T are the 
space-like and time-like extents of the loop, and the lattice vol- 
ume is L 4 (all dimensionful quantities are expressed in units of 
the lattice spacing a). A renormalized coupling can be defined 
by 



g 2 (R/L,L) = - 



R 2 



k(R/L) dRdT 



ln{W(R,T,L))\ T=R , (19) 



where for convenience the definition will be restricted to Wil- 
son loops with T = R, and (...) is the expectation value of some 
quantity over the full path integral. This definition can be moti- 
vated by perturbation theory, where the leading term is simply 
the bare coupling gl. The renormalization scheme is defined by 
holding R/L to some fixed value. The quantity k(R/L) is a ge- 
ometric factor which can be determined by calculating the Wil- 
son loop expectation values in lattice perturbation theory. The 
role of lattice simulations will be to measure non-perturbatively 
the expectation values. On the lattice, derivatives are replaced 
by finite differences, so the renormalized coupling is defined to 
be 



g 2 ((R + l/2)/L,L) = 
X(R+l/2,L) = -ln 



1 



-(R+l/2) 2 X (R + l/2,L), 



k(R/L) 

W{R + l,T + 1,L)W(R, T,L) 



\T=R : 



W{R+l,T,L)W(R,T + \,L) 

where x i s me Creutz ratio ||67l , and the renormalization 
scheme is defined by holding the value of r — (R + l/2)/L 
fixed. 

With this definition, the renormalized coupling g 2 is a func- 
tion of the lattice size L and the fixed value of r. The coupling is 
non-perturbatively defined, as the expectation values are calcu- 
lated via lattice simulations, which integrate over the full phase 
space of the theory. By measuring g 2 (r, L) non-perturbatively 
for fixed r and various L values, the running of the renormal- 
ized coupling is mapped out. In a QCD-like theory, g 2 increases 
with increasing L as we flow in the infrared direction. In a con- 
formal theory, g 2 flows towards some non-trivial infrared fixed 
point as L increases, whereas in a trivial theory, g 2 decreases 
with L. The advantage of this method is that no other energy 
scale is required to find the renormalization group flow. The 
renormalized coupling g 2 is also a function of the bare cou- 
pling g 2 y which is related to the lattice spacing a. Keeping the 
lattice spacing fixed, the running of g 2 (r, L) is affected by the 
lattice cut-off. The running has to be calculated in the con- 
tinuum limit, extrapolating to zero lattice spacing. A similar 
method was developed independently in [68]. 

One way to measure the running of the renormalized cou- 
pling in the continuum limit is via step-scaling. The bare lattice 
coupling is defined in the usual way (3 = 6/g^ as it appears 
in the lattice action. Some initial value of g 2 is picked from 
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: 1.636(23) 
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Figure 7: The method and the main test result for pure-gauge theory are 
shown in the figure. In the upper figure the extrapolation procedure 
picks up the leading a 2 /L 2 cutoff correction term in the step function. 
It gives the fit to the continuum limit value of the step function. In 
the lower figure, the running coupling g 2 (L) is shown. The blue points 
are from results on Creutz ratios using analytic/numeric Wilson loop 
lattice calculations in finite volumes with fixed value of r. In this pro- 
cedure we start from the one-loop expansion of Wilson loops in finite 
volumes based on the bare coupling |69|. The series is re-expanded 
in the boosted coupling constant at the relevant scale of the the Creutz 
ratio 1 70 1 to obtain realistic estimates of our running coupling with- 
out direct simulations. The rest of the procedure for the blue points 
follows what we described in the text. The green points are direct sim- 
ulation results, following our procedure. The running starts at the point 
g 2 (Lo) = 0.825. For almost all couplings there is excellent agreement 
with continuum 2-loop running. At the strongest coupling, the simula- 
tion results begin to break away from perturbation theory. 

which the renormalization group flow is started. On a sequence 
of lattice sizes L\,Li, ...,L n , the bare coupling is tuned on each 
lattice so that exactly the same value g 2 (r, Li, ft,) is measured 
from simulations. Now a new set of simulations is performed, 
on a sequence of lattice sizes 2L\,2Lq, 2L„, using the corre- 
sponding tuned couplings (3\,P2, ■■■,Pn- From the simulations, 
one measures g 2 (r, 2L/,/?,), which will vary with the bare cou- 
pling viz. the lattice spacing. These data can be extrapolated to 
the continuum as a function of 1/L,. This gives one blocking 
step L — > 2L in the continuum renormalization group flow. The 
whole procedure is then iterated. The chain of measurements 
gives the flow g 2 (r,L) -> g 2 (r,2L) -> g 2 (r,4L) g 2 (r,8L) 

as far as is feasible (Fig. [7]). One is free to choose a different 
blocking factor, say L — > (3/2)L, in which case more blocking 
steps are required to cover the same energy range. 

We applied the above procedure to the running coupling in- 
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P 


j 

L 


fermion mass 


trajectories 


8 (A) 


5 


12 


0.01 


318 


2.06(2) 




16 


0.01 


74 


1.67(11) 


7 


12 


0.01 


317 


1.207(5) 




12 


0.001 


116 


l. 207(12) 




16 


0.01 


198 


1.13(1) 


19 




01 


16? 






16 


0.01 


69 


0.577(9) 


15 


12 


0.01 


144 


0.447(3) 




12 


0.001 


91 


0.460(5) 




16 


0.01 


62 


0.444(7) 


25 


12 


0.01 


190 


0.255(1) 




16 


0.01 


156 


0.253(2) 



Table 1 : Running couplings bracketing the conformal fixed point of the 
Nf = 16 model in the conformal window. 

side the conformal window with Nf = 16 flavors. The shortcut 
of this pilot study ignores the extrapolation to the continuum 
limit. The running coupling therefore is still contaminated with 
finite cutoff effects. If the linear lattice size L is large enough, 
the trend from the volume dependence of g 2 (L, a 2 ) should indi- 
cate the location of the fixed point. For g 2 (L, a 2 ) > g* 2 we ex- 
pect the decrease of the running coupling as L grows although 
the cutoff of the flow cannot be removed above the fixed point. 
Below the fixed point with g 2 (L,a 2 ) < g* 2 we expect the run- 
ning coupling to grow as L increases and the continuum limit of 
the flow could be determined. The first results are summarized 
in Table 1. They are consistent with the presented picture. For 
example, at bare couplings [} = 5,7, 12 the cutoff dependent 
renormalized coupling is larger than 0.5 and decreasing with 
growing L. At small bare couplings the renormalized coupling 
is flat within errors and the flow direction is not determined. 
The independence of the results from the small quark mass of 
the simulations is tested in two runs at m q = 0.001. Precise 
determination of the conformal fixed point in the contiuum re- 
quires further studies. 
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